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For a wide class of technologically relevant compound III-V and II-VI semiconductor materials AC 
and BC mixed crystals (alloys) of the type A^Bi-zC can be realized. As the electronic properties like 
the bulk band gap vary continuously with x, any band gap in between that of the pure AC and BC 
systems can be obtained by choosing the appropriate concentration x, granted that the respective 
ratio is miscible and thermodynamically stable. In most cases the band gap does not vary linearly 
with x, but a pronounced bowing behavior as a function of the concentration is observed. In this 
paper we show that the electronic properties of such A^Bi-^C semiconductors and, in particular, 
the band gap bowing can well be described and understood starting from empirical tight binding 
models for the pure AC and BC systems. The electronic properties of the A^Bi^C system can be 
described by choosing the tight-binding parameters of the AC or BC system with probabilities x 
and 1 — x, respectively. We demonstrate this by exact diagonalization of finite but large supercells 
and by means of calculations within the established coherent potential approximation (CPA). We 
apply this treatment to the II-VI system Cd^Zni-^Se, to the III-V system In^Gai-^As and to the 
Ill-nitride system Ga^Ali-^N. 
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I. INTRODUCTION 

Disordered semiconductor alloys of the type A x Bi_ x C 
can be experimentally realized for many technologically 
important semiconductor materials AC and BC (like 
GaAs and InAs, CdSe and ZnSe, GaN and A1N and many 
more) by substituting the cation A by the B counter- 
part. For a wide class of materials, the value x can be 
adjusted over a large part of or even the whole concentra- 
tion range, depending on the miscibility gap of the alloy 
system. As the electronic and optical properties vary con- 
tinuously with x, the electronic properties (e.g. band gap 
and thus emission and absorption frequency, dielectric 
constant, etc.) can be tailored by choosing the appropri- 
ate concentration. Therefore, these semiconductor alloys 
have many applications, not only as bulk semiconduc- 
tors but also in low-dimensional structures like quantum 
wells, quantum wires and quantum dots. 

In this paper, we show that the electronic properties 
of such substitutional semiconductor alloys of the type 
A^Bi-^C can be described and understood on the basis 
of empirical tight-binding models (ETBM) for the pure 
semiconductors AC and BC. Once the tight-binding (TB) 
parameters of AC and BC are known, the ETBM for the 
substitutionally disordered Aa,Bi_ x C alloy is obtained by 
choosing the AC TB parameters with probability x and 
the BC TB parameters with probability 1 — x. The elec- 
tronic eigenenergies and thus the density of states and 
the band gap of this disordered system can then be deter- 
mined by exact diagonalization of a finite (but large) su- 
percell. Unfortunately, it is necessary to consider a large 
number N of microscopically distinct configurations for 
each concentration well as sufficiently large super- 

cells, consisting of several thousand primitive unit cells 



at least. As this procedure is numerically expensive, we 
also considered and applied well established approxima- 
tions for the treatment of disordered systems, namely the 
coherent potential approximation (CPA) and the virtual 
crystal approximation (VCA). While the CPA is a self- 
consistent Green function approximation, the much sim- 
pler (but still frequently applied) VCA corresponds to a 
static mean-field treatment of the disorder by replacing 
the randomly fluctuating potential by an average poten- 
tial. 

We apply our treatment to the zincblende phase of 
three interesting materials, namely to the II-VI semicon- 
ductor alloy system Cd^Zni-^Se, to the technically es- 
tablished III-V alloy In^Gai-^As and to the Ill-nitride 
system Ga^Ali^^N. In a previous paper p] it has al- 
ready been shown that the supercell exact diagonaliza- 
tion method yields excellent agreement with experimen- 
tal results obtained for Cd^Zni-^Se. But as [T] was a 
joint paper together with experimentalists, the theoreti- 
cal method was not yet described and discussed in detail 
and no comparison with CPA and VCA results for this 
material has been presented there. Whereas the VCA, 
in general, yields no band gap bowing, the CPA may 
even overestimate the bowing, as we will show in the 
present paper. We particularly present results obtained 
within combinations of CPA and VCA, treating the site- 
diagonal disorder in CPA and the off-diagonal disorder 
in VCA. Furthermore, we explicitly discuss the influence 
of a properly chosen basis set for the bowing results. 

Also for In^Gai-^As, reasonable results for the band 
gap bowing are obtained in agreement with literature val- 
ues as obtained by experiment. For Ca^Ali^N the band 
gap bowing is rather small and the VCA results are not 
too bad at all; a crossover from a direct band gap to an 
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indirect band gap behavior as a function of x is obtained 
in all three applied approaches. 

This paper is organized as follows: In Section [TTJ we 
present the necessary theoretical background for our cal- 
culations of the electronic properties of A^Bi-^C semi- 
conductor bulk alloys. It contains a brief description of 
the employed tight-binding model, the CPA and the mod- 
elling on a finite ensemble of supercells. Section |III| con- 
tains a detailed analysis of the results that we obtained 
for the electronic properties of Cdo.5Zno.5Se with the dif- 
ferent models. Furthermore, we carefully analyze the re- 
sulting band gap bowing of Cd a; Zn 1 _ ;r Se, In^Gai-^As 
and Ga x Ali_ x N and discuss the applicability of the CPA 
and VCA to these systems. 



II. THEORY 



Empirical tight-binding model for pure bulk 
semiconductors: 



We know from textbook solid state physics that one 
complete basis set of the translationally invariant pure 
bulk system is given by the Wannier states |nR), where 
n is the band index and R denotes the unit cell where 
the state is predominantly localized. The Wannier states 
are connected to the itinerant Bloch states |nk), where 
k is any wave vector in the first Brillouin zone (BZ), by 
means of a discrete Fourier transformation. The ETBM, 
as originally introduced by Slater and Koster [5J, uses a 
finite linear combination of atomic orbitals on the basis 
sites as ansatz for the Wannier states. The matrix ele- 
ments are then fitted in order to sufficiently reproduce 
prominent band structure features. 

Because the atomic states, and hence the Wannier 
states, are not explicitly used in an ETBM, one has the 
additional freedom to directly assume a finite basis set 
of localized Wannier-like states. As the Wannier basis 
carries the translational invariance of the crystal, these 
states merely have to be assigned to the sites R of the 
underlying Bravais lattice. This parametrization scheme 
is also known as effective bond-orbital model (EBOM) in 
the literature. 

In this paper, we will use a sp 3 basis per spin direc- 
tion {t, 4} to reproduce N a — 8 bands for crystals with 
zincblende structure: One spin-degenerate s-like conduc- 
tion band (CB) and three spin-degenerate p-like valence 
bands (VB), in detail the heavy hole (HH), light hole 
(LH) and split-off (SO) band. Thus, R labels the N R 
sites of the underlying fee Bravais lattice, and the TB 



matrix elements of the bulk Hamiltonian H 

E™' = (Ra|ff bulk |RV), 
with the orbitals 



bulk 



(1) 



The translational invariance in the pure crystal reduces 
the N a Nn, x N a Nn matrix problem to a N a x N a problem 
for each k [2J, which will not be the case for disordered 
systems and in the later introduced supercell method. 
The band structure E(k) is then given by the solution of 
the eigenproblem 
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e lk - R ^ CQ '(k) = £;(k) CQ (k). 
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If one restricts the non-vanishing matrix elements, 
Eq. up to a finite neighborhood and makes proper 
use of symmetry relations, the TB matrix elements can 
be expressed as a function of a set of material param- 
eters of either the AC or BC material (e. g. CdSe and 
ZnSe). We will use the parametrization scheme of Loehr 
[3], which includes coupling up to second nearest neigh- 
bors. The parameters are given in the appendix. Figure 
[T] shows the band structure and density of states (DOS) 
as calculated in the sp 3 EBOM for cubic CdSe and ZnSe, 
respectively. 

Despite the small basis set, this TB model notably 
allows for the reproduction of a realistic bandstructure 
and bandwidth throughout the whole BZ. This is a typ- 
ical feature of parametrizations with Wannier-like basis 
states [THj and makes the EBOM especially suitable for 
the purpose of the present work, as both the CPA and 
the supercell calculations for the A^Bi^C alloys will re- 
quire a realistic input DOS for the pure cases x = and 
x = 1 as starting point — the relative position of the band 
centers and the bandwidths will crucially influence the 
alloy properties. We have also already mentioned in the 
introduction that the inclusion of further bands would 
significantly increase the computational effort. Finally, 
the resolution on the scale of unit cells instead of atomic 
sites will turn out as favorable for the simulation of the 
ternary A^Bi-^C materials, as it will allow for an unam- 
biguous assignment of the alloy lattice sites to either AC 
or BC (see Sec. |IlD]). 
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B. The coherent potential approximation 

The basic principle of the CPA has independently been 
developed by several groups in the late sixties of the last 
century 4] [5 [6 . The best known formulation is certainly 
the work of Soven 0], which explicitely dealt with the 
calculation of the electronic DOS of substitutionally dis- 
ordered one-dimensional systems. An excellect and com- 
prehensive general introduction into the CPA formalism 
can be found in [7J. 

The simplest form of the CPA for an alloy assumes 
uncorrelated substitutional disorder of the respective 
species. If all sites are indistinguishable (i. e. no divi- 
sion into sublattices is necessary), as it is the case in the 
spatial discretization in the EBOM, the probability of 
finding either species is just given by the concentrations 
x and 1 — x. 
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FIG. 1: (Color online) Band structure (left) and DOS (right) of CdSe and ZnSe, calculated in the sp 3 EBOM. The energies of 
CdSe are already shifted by a valence band offset A_E v b = 0.22 eV with respect to ZnSe. For later use, the energetic position 
of the first moments of the bands is also given. See text and appendix for details. 



For the sake of clarity, we will supress the band and 
orbital indices until further notice. In order to be consis- 
tent with the common CPA literature, let us denote the 
site diagonal (R = 0) TB matrix elements of the AC or 
BC material as t>Ac/BC : 



V AC/BC 



E' 



00 

AC/BC 



rrbulk 

-"AC/BC 



We will now briefly specify the most important assump- 
tions that enter the CPA (see e. g. [7] for details of the 
derivation) : 

1. The disorder is confined to the diagonal elements. 
The off-diagonal (or hopping) matrix elements with 
R R' are either identical for both species or 
can approximately be replaced by a common value, 
e. g. the VCA average. If we now separate the TB 
Hamiltonian of the alloy H into a diagonal part V 
and an off-diagonal part W, such that H = V + W, 
only V = 5Zr,^ r |R)(R| is site-dependent, as 
v = wac/bCi depending on the species on the site 
R. The operator W is still translationally invari- 
ant under translations by R and therefore remains 
diagonal with respect to Bloch states |k). 

2. The configurational average (...) over the resolvent 
of H defines an effective Hamiltonian H e g(z): 



((zt-H) 
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\zt 
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Here, 1 is the identity operator and z is in the 
complex energy plane, containing the energy axis 
E = Re(z). The self-energy operator £(z) = 
Hefi(z) — W absorbs the influence of the disorder on 
the microscopic scale. Note that H e g will in general 
be non-hermitian. 



3. Due to the single-site nature of the CPA, the self- 
energy is diagonal in every representation. Further- 
more, its matrix elements S°(z) are neither depen- 
dent on k nor on R. 

In order to obtain the self energy matrix elements 
which uniquely define the effective medium, we have to 
solve the CPA equation for two constituents: 



X \ V K Q 



(1 - x) [v BC - YP{z)} 



1 - [v AC - E°(*)] G R (z) 1 - [v BC - E°(«)] G R ( 2 



= 0. 

The complex- valued Gr(z) is the configurationally aver- 
aged one-particle Green function in Wannier representa- 
tion (more precisely its R-diagonal element). Although 
Gr(z) does not depend on R in the CPA, we will keep 
the index to clearly distinguish it from its Bloch repre- 
sentation G-k(z). The two representations are connected 
by 



G R (z) = 



EqJZ} 
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with Ak as the number of wave vectors in the first BZ 
or the corresponding irreducible wedge (as Ak = Ar 
and Gr = Gr+r' , this relation follows directly from the 
invariance of the trace under unitary transformations). 

Under certain conditions, the summation over all k- 
valucs can be avoided by the introduction of a sufficiently 
smooth DOS gw for the off-diagonal part W: 

9w(E) = ~J2 S ( E -W W \Vy (7) 

k k 
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Now the calculation of Gr(z) can be performed as a one- 
dimensional integration/summation over the energy axis: 



G R (z) 
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]9w{Ei) J dE 
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£,: 
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In [z - E°(z) - E l+1 ] } 



The approximation holds if gw (E) is approximately con- 
stant on the interval [Ei, E^±]. 

The equations ^ and optionally Q or pi) can be 
solved in a self-consisted manner in order to obtain E°(z) 
and Gr(z). We can then directly obtain further quan- 
tities for the effective medium; e. g. the configurationally 
averaged DOS g(E) per lattice site is given as 



9(E) 



— lim Im Tr 

7T z-*E+ N R 

1 



{[zl-Htf(z)] 



-- lim ImG R (£; + iS). 

7T S\0 
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As a rule of thumb, the CPA is known to yield very 
good results in two limit cases [7J: 

1. The weak scattering limit, where the difference of 
the first moments of the substituents (aka the cen- 
ters of gravity of the bands) is smaller than the 
respective bandwidths. 

2. The split band limit or atomic limit, where the sub- 
stituents' first moments are either sufficiently far 
apart or the bandwidths are small, so that the re- 
spective bands do not overlap. 

Additionally, it of course also gives the correct pure limit 
for x — > and x —> 1 . Also, in the limit case of vanishing 
difference of the moments, the CPA reduces to the VCA. 



C. Multiband CPA + ETBM 

The formal extension of the CPA to multiband TB 
models is straightforward and has been used in differ- 
ent levels of detail to qualitatively examine the elec- 
tronic properties of disordered alloys (see e. g. [H][5] for 
Si.Gei-,, [HI for Pc^Hi-*, (Tp for Cd.Hg^Te or QJ] 
for palladium-noble-metal alloys). 



1. sp representation 

In the sp 3 EBOM, the localized basis is now given by 
the TB orbitals |Ra). Hence, Eq. ([5| has to be replaced 



by the corresponding matrix equation, with 



^ac/bc - 
E°(js) -> Y, sp \z 



AC/BC 
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G R (z) -> G% (z) = 



aa' AC/BC 



G R Q (z) 



(10) 



as 4 x 4 matrices per spin direction (here, the square 
brackets denote matrices). The effective Hamiltonian 

matrix (z) in the sp 3 TB scheme can be obtained 
by the substitutions 



CPA 



(*), 



pRR' VC A pRR' 



(11) 
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sJ^bc (12) 



i. e. the hopping matrix elements are approximated in 
the VCA. Like in Eq. ^ , we then obtain the Green func- 
tion via BZ summation, 



HS (*) 



(13) 



where 1 is the 8x8 identity matrix. The matrix G R (z) 
has non- vanishing off-diagonal elements G R Q (z), as the 
TB Hamiltonian is not diagonal in the orbital basis |Ra). 



2. Band-diagonal (Wannier) representation 

Alternatively, we can again calculate Green's function 
via energy integration. In a band-diagonal (Wannier-like) 
TB representation, the site-diagonal TB matrix element 
E° = (On |iJ bulk | On) of the n-th band can be calculated 
from the first moment 



dEEg n {E), 



(14) 



with g n (E) as the DOS of band n (normalized to unity). 
Therefore, we can assign the moments of the AC or BC 
system to a corresponding self-energy: 



E 



: AC/BC 



TT{z). 



(15) 



The analogon of equation ^ is then the band-diagonal 
Green function matrix G R n (z) (denoted by "wn" for 
"Wannier" from now on) with elements 

Gr'W « Y,9w(E n ^){\n[z-^{z)-E n ^ 

i 

-In [z-Y. n {z)-E n , i+1 ]}5 nn ,{l&) 

with gfo(E) = g VCA (E - E°) as the VCA DOS of the 
ri-th band, shifted by E®. This procedure allows us to 
calculate the DOS of the off-diagonal operator W in the 
Wannier basis by diagonalization of the EBOM Hamilto- 
nian in the VCA. Notably, we do not have to explicitly 
know the corresponding matrix elements (Rn \ W\ R'n). 
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3. Calculation of electronic properties 

Finally, the CPA + ETBM DOS of the alloy can be 
calculated by either tracing over the orbital index (in the 
sp 3 basis) or the band index (in the Wannier basis) of 
the matrix elements: 

9 sp3 (E) = -^limlmYsG^iE + iS), (17) 

a. 

9 wn (E) = -^limlmJ2G^(E + i6). (18) 

n 

Although the trace of a matrix is of course independent 
of the representation used, it should be noted that g sp 
and g wn are not identical: The self-energies, which de- 
fine the effective medium, replace different quantities [see 
Eqs. (11) and (15)], which results in different operators 
H eS . 

Due to the translational invariance of the self energy 
matrix elements and thus H c g, we can define a complex 
band structure of the medium. The electronic excitations 
of the effective medium can be assigned to quasiparticles 
with a modified dispersion relation 



Ej?(k) = (kn \W\ kn) + Re E n (£^). 



(19) 



The corresponding CPA one-particle spectral function 
1 



luaImGZ n (E + i5), (20) 

7T S\0 



will in general be broadened in case of disorder (0 < x < 
1), as the non- vanishing imaginary part of S n accounts 
for a finite lifetime r cx 1/ImE™ (E) of the excitation. 

As an example, Fig. [2] visualizes the complex band 
structure E^(k) of Cdo.5Zno.5Se on the left, calculated 
in the combination of the CPA and the sp 3 EBOM. On 
this energy scale over the whole bandwidth, the results 
from the Wannier and the sp 3 representation are not 
distinguishable. The color-coding clearly shows the In- 
dependent broadening of the band structure due to fi- 
nite lifetime effects. The figure on the right additionally 
shows the corresponding CPA DOS and the spectral func- 
tion at the BZ center for all four spin-degenerate bands, 
i. e. 



S^o(E) 



1 



7T S\0 



Note that we have chosen a relatively large imaginary 
part S = 0.05 eV for the energy axis, which leads to a 
small but non-vanishing DOS in the band gap region. 
We furthermore used #k ss 10 6 values in the irreducible 
BZ of the fee lattice. The peak structure of the k = 
quasiparticle excitations at the edges of the band gap 
region is clearly visible. A qualitative and quantitative 
discussion of the electronic structure in comparison with 
VCA and supercell results will follow in Section III A 



D. Supercell tight-binding calculation in the 
EBOM 



If the Hamilton operator H is no longer translation- 
ally invariant, the TB approach is not reducable to the 
form (|3| and we are left with the N a N-R x N a N-R matrix 
eigenvalue equation 



£ (Ra|#|RV)cv R , =Ec, 



(21) 



'R' 



We will now use a finite supercell with periodic bound- 
ary conditions. Like in calculations for zero-dimensional 
nanostructures [13] [14] [15] [16] , we will furthermore as- 
sume that 



(Ra \H\ RV) « (Ra |ff bulk | R'a') = E, 



tiRR' 



(22) 



Again assuming uncorrelated substitutional disorder in 
the A^Bi-zC alloy, each primitive cell will be occupied 
by the AC or BC basis, where the probability of finding 
either an AC or BC pair in the unit cell at R is directly 
given by the concentrations x and 1 — x. Building the 
Hamilton matrix for the supercell, we will therefore use 
the matrix elements of the pure AC or BC material for 
the corresponding lattice sites. The valence band off- 
set between the two materials is incorporated by shifting 
the respective site-diagonal matrix elements by a value 
A£ v b. Hopping matrix elements between unit cells of 
different material are approximated by the arithmetic av- 
erage of the corresponding AC or BC values. Although 
this approach for the incorporation of the energy offset 
and the hopping between two materials is very simple, 
it gives the correct limit in the pure case. In the phase 
separation case (a limit of which the propability is prac- 
tically zero in case of uncorrelated disorder), it would 
furthermore lead to an interface treatment that has been 
extensively tested in calculations for low-dimensional het- 
erostructures [H] [H] [IS] . Specifically, the results turn 
out to be insensitive to small variations of the hopping 
between two materials (e. g. the usage of a geometric 
instead of an arithmetic average). 

The usage of the EBOM, i. e. the usage of Wannier-like 
orbitals situated at the sites of the Bravais lattice in the 
supercell approach, has several advantages over similar 
approaches which use an ETBM with discretization on 
atomic sites: 

1. Even when introducing disorder on a microscopic 
scale, each lattice site can be unambigously as- 
signed to one site-diagonal TB matrix element and 
the corresponding band structure parametrization 
for either AC or BC. In ETBM supercell calcula- 
tions with atomic resolution, each anion of the type 
C will locally be surrounded by a different number 
of A or B cations, thus making an assignment of 
the diagonal elements of the anions to the band 
structure of either AC or BC impossible. 

It is common to then use either a concentrationally 
averaged VCA value or to determine this matrix 



6 




FIG. 2: (Color online) Complex band structure (left), DOS and spectral function at k = (right) for Cdo.5Zno.5Se, calculated 
in the combination of the CPA and the sp 3 EBOM. To enhance contrast, we have color-coded the logarithm of — ImGk- The 
spectral function and the DOS are normalized to a common scale. 



element as a weighted average of the C matrix ele- 
ments for AC and BC, depending upon the number 
of nearest-neighbour atoms A or B [T7][IH]- This 
will lead to an effectively more coarse-grained res- 
olution as in the case of the EBOM, as the latter 
model only has to average the intersite hopping ma- 
trix elements (which typically differ on a scale of 10 
meV in materials with moderate lattice mismatch) . 
In the Cd^Zni-^Se system for example, the site di- 
agonal matrix elements for the Se anions in CdSe 
and ZnSe (see e. g. Refs. [T7] and [T!5]) differ to a 
larger extent than the EBOM hopping matrix ele- 
ments when using a congruent set of input param- 
eters. 

2. All influences which originate from effects on a 
smaller length scale, like the difference in the AC 
and BC bond lengths, are absorbed into the values 
of the corresponding TB matrix elements between 
the effective orbitals. 

3. As the results are not sensitive to the exact treat- 
ment of the hopping between AC and BC sites, fur- 
ther effects that basically result in minor variations 
of the hopping matrix elements (like small bond an- 
gle changes due to relaxation) can be neglected in 
a first approximation. 

The numerical diagonalization (e. g. using standard nu- 
merical libraries like ARPACK/PARPACK) of the corre- 
sponding Hamiltonian for a fixed concentration and a fi- 
nite number of microscopically distinct configurations 
gives the density of states (DOS) of the finite ensemble. 
In order to obtain a meaningful DOS from the super- 
cell calculations, we must first appropriately define it. 
Strictly speaking, a macroscopic alloy crystal represents 



just one realization; by dividing it into small portions, 
we can nevertheless get subsystems that differ from each 
other on a microscopic scale. In this sense, the DOS for 
one fixed concentration x is then obtained by the average 
of the DOS for each finite ensemble. 

To eliminate the influence of finite size effects, the num- 
ber of lattice sites Nn as well as the ensemble size N must 
be sufficiently large. We point to previous work [1] for a 
careful analysis of the convergency behaviour and will use 
the recommendations throughout this paper. In a nut- 
shell, a resolution of the band edges up to 0.01 eV, which 
is the typical input accuracy for the material parameters, 
will require supercells with ATr w 2000-4000 lattice sites 
and N fa 50 microscopically distinct configurations per 
concentration. For a discussion of properties on a larger 
energy scale, e. g. a comparison of the DOS over the 
whole bandwidth, smaller supercells and ensemble sizes 
can be chosen. 

In contrast to the CPA, the supercell approach can 
easily be augmented to simulate effects not only of con- 
figurational, but also of concentrational disorder. This 
can easily be achieved when we drop the constraint that 
the overall concentration of AC sites N^ c /N^ per config- 
uration i should equal the point probability x and occupy 
each lattice site independently. In the limit of large N, we 

will of course have limjv-j-oo S,=i ^Yac/^Yr = x -> so that 
the constraint is fulfilled for sufficiently large ensemble 
numbers. For most cases, this is closer to experimental 
reality anyway, as concentration values are commonly av- 
erages over macroscopic volumes, e. g. by means of X-ray 
diffraction pQ. It also allows us to perform calculations 
for concentration values x where N-r/x (fc IN. 
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FIG. 3: (Color online) Direct comparison of the DOS of 
Cdo.5Zno.5Se, calculated in the combination of the VCA 
(blue), CPA (red) and supercell approach (black) and the 
EBOM. See text for further details. 



III. RESULTS 

In this section, we will apply the CPA EBOM and the 
supercell EBOM to cubic Cd^Zni^^Se, In^Gai-^As and 
Ga-jAli-^N. In addition, we will also add results that are 
obtained by a pure VCA calculation. 

Besides the fact that we have already shown the reli- 
ability of the supercell EBOM for cubic Cd^Zni-^Se in 
comparison with experimental results in [1] (albeit for 
slightly different material parameters to meet the exper- 
imental boundary conditions), this material systems is 
also especially interesting for a quantitative and quali- 
tative analysis of the applicability of the CPA. A closer 
look at the DOS and the energetic position of the first 
moments in Fig.[T] reveals that the conduction bands of 
CdSe and ZnSe neither fulfill the weak scattering nor the 
split band condition very well, as the energetic range of 
the overlap is comparable to the difference of the first mo- 
ments. The In^Gai-^As system under consideration will 
in contrast be closer to the weak scattering limit. This 
condition will also apply to the zincblende Ga^A^^N al- 
loy, which additionally comprises a direct-indirect band 
gap transition at a certain mixing ratio. 



Overall, the alloy DOS of the valence bands is very 
similar in all three models. This is not very surprising, 
as the substitutional disorder is restricted to the cations 
of the material, and the valence bands mainly stem from 
atomic p-orbitals of the Se anions. However, the conduc- 
tion band DOS accordingly shows different features in 
the three models. 

It is clearly visible that the VCA DOS is an interpo- 
lation of the DOS of the pure CdSe and ZnSe material 
as depicted in Fig. [I] aside from an energetic shift of the 
bands, no new features arise for the alloy material. 

Contrary to the VCA, the CPA gives a DOS with 
qualitative and quantitative features that exceed the re- 
sults of simple interpolation schemes by far. The con- 
duction band visibly splits into two subbands. When 
the resolution is further increased, a quasi-gap can be 
identified. The relative spectral weight of the two sub- 
bands is exactly given by the concentration ratio of 
x/(l — x) = 0.5/0.5 = I of the substituents (this also 
holds for all other values of x with identifiable subband 
splittings). 

Like in the CPA, the supercell EBOM also gives an al- 
loy DOS of which the structure is more complicated than 
the DOS of the constituents. The conduction band DOS 
again splits into a two-subband structure, where the rela- 
tive spectral weight equals the concentration ratio of the 
substituents. Furthermore, the two subbands can now 
clearly be distinguished by a difference in their shape. 
Additionally, the quasi-gap is located at a lower energy 
than in the CPA. 

Overall, the fact that VCA totally fails to reproduce 
the additional features of the alloy's conduction band was 
to be expected; the one-electron potential which enters 
the hopping matrix elements is not a self-averaging quan- 
tity (in the sense that it can be replaced by its ensemble 
average for a sufficiently large sample) , as opposed to the 
one-electron Green's function [20] , 

The artificial symmetry in the conduction subband 
structure in the CPA is an artefact of the usage of concen- 
trationally dependent, but nevertheless common intersite 
hoppings. This mean- field approach for the hopping ma- 
trix elements directly carries over to the shape of the DOS 
(this can most easily be seen in the Wannier basis, as the 
shape and the bandwidth is eventually determined by the 
hopping elements, while the band moments are given by 
the site-diagonal elements, see Sec. II C 2). 



A. Comparison of the overall density of states of 
Cdo.5Zno.5Se 

In this section, we will compare the overall DOS of the 
A^Bi-a-C alloy as calculated with the CPA and the su- 
percell EBOM, using the example of Cdo.5Zno.5Se, along 
with results from the simple VCA. The supercell calcu- 
lations were performed with 20 microscopically distinct 
configurations and on cubic supercells with 2048 lattice 
sites, i. e. 4096 atoms; the numerical parameters for the 
CPA calculation match those of Fig. [2j 



B. Comparison of the band gap bowing of 

Cd^Zm-^Se 

We will now extensively examine the accuracy of the 
CPA approach for properties on a smaller energy scale 
and use the example of the single-particle band gap for 
a qualitative as well as quantitative analysis. 

Most alloyed bulk semiconductors show a more or less 
pronounced bowing of the band gap E g as a function 
of the concentration x. The simplest way to describe 
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FIG. 4: (Color online) Conduction band edge of Cd^Zni-^Se, 
calculated in the CPA. The peak of the spectral function at 
the BZ center coincides with the band edge within the energy 
resolution. See text for further details. 



the deviation from a linear behaviour is the assumption 
of a parabolic E g (x) curve and therefore the use of a 
single, concentration-independent bowing parameter b, 
such that 

E s (x) = xE^ c + (1 - x) E* c - x(l - x)b. (23) 

Here, the indices AC and BC assign the properties of 
the pure binary materials. In general, the literature val- 
ues for b show a surprisingly large variety even for ap- 
parently comparable experimental conditions (the reader 
may check comprehensive review articles like [21] or |22j). 
For the II- VI bulk alloy Cd :! .Zn 1 _ x Se for example, a broad 
range of values between 6 = and b = 1.26 eV has 
been reported throughout publications from the last two 
decades [H] [23] [24] [25] . The large disparity on the 
experimental side can for example result from difficult 
growth conditions for the mixed systems. On the theo- 
retical side, the inadequate use of too simple approaches 
like the VCA can lead to wrong results. 

In the CPA, the band gap can in principle be read off 
from the DOS. For numerical reasons, the CPA DOS will 
not completely fall to zero in the bandgap, as a finite 
imaginary part 8 of the energy is required. Nevertheless, 
it is possible to identify the band gap with desired accu- 
racy by increasing the resolution. We used an imaginary 
part of 6 = 1CP 4 eV and #k ss 10 7 values in the irre- 
ducible BZ. The high BZ resolution has turned out to 
be crucial to obtain convergence for the results for the 
band gap. In the Wannier representation, it is further- 
more very important to carefully discretize the VCA DOS 
when using Eq. ( 16 ) because the DOS will contain kinks 



that stem from Van Hove singularities [critical points 
where grad k _E„(k) = 0]. The usage of TB models that 
give a reliable band structure throughout the whole BZ 
(in contrast to dispersions from effective mass or k • p 
models) will additionally lead to sharp peaks at some of 



these critical points, as the slope of non-degenerate bands 
must also vanish at the BZ boundaries. 

Figure [4] shows the conduction band edge region of 
Cdo.5Zno.5Se, calculated in the Wannier CPA. The peak 
of spectral function of the k = quasiparticle obviously 
coincides with the band edge within the chosen energy 
resolution of Az = 5 x 1CP 3 eV. A further analysis 
(not shown) reveals that the CPA band edge states are 
mostly (5-like; in case of broadened peaks, the correspond- 
ing linewidth is so small that the peak position still allows 
for a convenient determination of the conduction band 
edge E c and the valence band edge E v with an accuracy 
of 0.01 eV. Hence, the band gap E g (x) = E c (x) — E v (x) 
of the alloy can be determined with an accuracy of 0.02 
eV. 

The corresponding supercell band gap is given by 



E s (x) = min{Ei(x)} 



:{El(x)}, (24) 



where i again numbers the distinct configurations. Note 
that this definition of the band gap implies 



E g (x) + 



-y. 



[Ei(x) - El(x)] , 



(25) 



i. e. the band gap of the disordered alloy is different from 
the configurational average over the energy gaps for fixed 
realizations i. 

The resulting curves for the Cd^Zni-^Se band gap are 
depicted in Fig. [5] In order to perform a detailed exam- 
ination of the influence of the disorder and the applica- 
bility of the CPA, we performed two different supercell 
calculations for the Cd^Zni-^Se alloy system. In order 
to get rid of finite size effects, we used supercells with 
Vr = 4000 lattice sites and calculated the band edges for 
N = 50 distinct configurations per concentration. The 
concentration itself is varied in steps of 0.1. 

The left subfigure contains results from the supercell 
EBOM exactly as described in Sec. |IID[ thus including 
the full disorder in the site-diagonal and the hopping ma- 
trix elements on the microscopic scale. The supercell 
results in the right subfigure have been obtained under 
the artificial restriction to site-diagonal disorder. This 
means that only the site-diagonal matrix elements differ 
throughout the cell and with each configuration, while 
the hopping matrix elements for each concentration were 
substituted by their VCA values. Both the supercell ap- 
proach and the CPA then use the same level of mean-field 
approximation for the hopping matrix elements. As each 
of supercell bowing curve requires the partial diagonal- 
ization of 9 x 50 = 450 Hamiltonian matrices (the x = 
and x = 1 values are input material parameters), we ad- 
ditionally neglected the spin-orbit coupling for this com- 
parison. Besides an increase in computation time by a 
factor 2 3 = 8, the memory comsumption is lowered, as all 
matrix entries are then real numbers. We also added the 
pure VCA results to both figures — as the tight-binding 
matrix is diagonal at k = 0, we are left with a linear 
E g (x) curve, which is at least able to indicate the devia- 
tion of the other results from a linear interpolation. 
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FIG. 5: (Color online) Concentration-dependent band gap E g (x) of Cd^Zni-^Se, calculated in the VCA (blue), sp 3 CPA 
(magenta), Wannier CPA (red, only in left figure) and supercell method (black). The spin-orbit coupling is set to zero. In the 
left figure, the supercell results map the full disorder. In the right figure, the disorder in the supercells is artificially restricted 
to the site-diagonal elements in the sp 3 basis. See text for more details. 



We will first turn to the full disorder case on the left. 
Overall, the bowing obtained in the CPA is clearly larger 
than the corresponding supercell result. When the Wan- 
nier representation is used, the bowing curve is slightly 
closer to the supercell case than in the sp 3 basis. This 
was to be expected, as the overall shapes and bandwidths 
of the pure CdSe and ZnSe bands are very much alike 
(see again Fig.[lJ. Nevertheless, the CPA (as well as the 
VCA) fails to satisfactorily reproduce the supercell band 
gaps over the whole concentration range. We can even 
identify a slight "overbowing" , where the CPA band gap 
of the alloy dips beneath the value for pure CdSe for the 
Cd-rich concentrations. This effect is ultimately a conse- 
quence of the usage of the concentration-dependent VCA 
average for the non-diagonal part of the Hamilton. Con- 
sequently, it does not occur in standard textbook exam- 
ples, where the same hopping values for the constituents 
are used and one is left with a concentration-independent 
bandwidth. 

A look at the right subfigure, where the supercell re- 
sults with diagonal disorder are given along with the sp 3 
CPA curve, clearly reveals the importance of the non- 
diagonal part (as the supercell potential lacks transla- 
tional invariance, our supercell TB model can only be 
used in the sp 3 basis, so that the Wannier CPA results 
are also not given again here). We easily notice that the 
CPA and supercell results now coincide very well. Con- 
sequently, we can state that the CPA can in principle 
simulate the influence of the disorder in the site-diagonal 
elements on the band gap and the deviations stem from 
the mean-field treatment for the hopping. If we enforce 
the constraint N^/N^ = x for single concentration val- 
ues and thus eliminate the small influence of concentra- 
tional disorder (see Sec. II D) on the supercell results, 
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the discrepancy is even smaller (not shown). 

For the sake of completeness, we depict the Wannier 



FIG. 6: (Color online) Concentration-dependent band gap 
E s (x) of Cd^Zni-zSe including spin-orbit coupling, calcu- 
lated in the VCA (blue), Wannier CPA (red), and supercell 
method (black). 



CPA and supercell results including spin-orbit interac- 
tion in Fig. [6] (the sp 3 CPA results will be omitted from 
now on). We also added errorbars that account for the 
reading accuracy and finite size effects. The difference 
in the first moments of the CdSe and ZnSe conduction 
bands turns out to be slightly smaller when the spin is 
included. Obviously, the bowing is reduced and the de- 
viation between the Wannier CPA and supercell results 
decreases further. Still, only the x > 0.9 results over- 
lap within the error range. The best possible second- 
order fit to the E g (x) curve yields a bowing parameter 
of b = (0.71 ± 0.09) eV for the supercell results and 
b = (1.0 ± 0.1) eV for the curve as calculated with the 
Wannier CPA (note that the bowing values in Q] were 
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calculated for slightly differing band gap values in the 
pure case) . It should be emphasized that the error range 
for the bowing values can only account for the influence 
of the reading accuracy and finite size effects, as well as 
for the deviation from a parabolic behaviour (this is al- 
ways to be expected for differing lattice constants, see 
|26|). The disparities that can arise from the uncertainty 
in the input parameters for the pure materials (band 
gaps, effective masses, valence band offsets) from differ- 
ent sources in the literature cannot reliably be estimated 
at reasonable expense, so that the error range is certainly 
underestimated. 



C. Band gap bowing of In^Gai-^As 

As a second example for the calculation of the band 
gap bowing, we apply the methods to the zincblende bulk 
alloy In^Gai-zAs. The band structures, DOS and cor- 
responding first moments of InAs and GaAs are given in 
Fig.[7] A valence band offset of AE vh = 0.85 eV [27] has 
already been incorporated. This value has been obtained 
by the relative energetic position to transition-metal im- 
purities and is in agreement with experimental data from 
measurements on Au Schottky barriers |28) . 

For this parameter set of InAs and GaAs the first mo- 
ments of the conduction band are closer together than in 
the case of CdSe and ZnSe, going along with a larger over- 
lap of the bands. Hence, we are closer to the weak scatter- 
ing condition and therefore expect better results for the 
concentration-dependent band gap from the (Wannier) 
CPA. Again using N = 50 configurations and Ar = 4000 
lattice sites for the supercells and also the same numerical 
parameters for the Wannier CPA, we obtain the bowing 
curve depicted in Fig. [8] 

Overall, the deviation from a linear behaviour (again 
indicated by the VCA results) is smaller for In^Gai-^As 
than for Cd^Zni-^Se. Furthermore, all results from the 
supercell calculations and the CPA calculations now over- 
lap within the error range, although the CPA still slightly 
overestimates the bowing behaviour. A second-order fit 
of the E g (x) values gives a bowing of b = (0.39±0.04) eV 
for the supercell results and b = (0.49 ± 0.05) eV for the 
curve as calculated with the Wannier CPA. The fact that 
the resulting ranges touch but do not overlap is mainly 
due to the deviation from a parabolic curve. Neverthe- 
less, including the error boundaries both values lie in the 
range between b = 0.32-0.46 eV that is recommended 
in the literature by Vurgaftman et al. in |21j . More re- 
cent ab initio calculations within the DFT+LDA (which 
are known to systematically underestimate the band gap) 
give a slightly larger bowing in the range of 0.5-0.8 eV^H 

D. Band gap bowing of Ga^Ali-^N 

As a final example, we will calculate the concentra- 
tion dependent band gap of the zincblende phase of 



Ga^Ali-zN, which is frequently used as barrier mate- 
rial in optoelectronics [21] . While the wurtzite modifica- 
tion of A1N is the only Al-containing III-V semiconduc- 
tor with a direct band gap, its zincblende modification 
is most likely indirect with the conduction band mini- 
mum at the A-point and the valence band maximum at 
the T-point [301 (although a direct band gap is sometimes 
also assumed, see e. g. |31[) . As can be seen in Fig.[9) the 
sp 3 EBOM is able to reproduce the indirect band gap 
of A1N properly, as it is also fitted to the A-point en- 
ergies. We use a valence band offset of AE vh = 0.8 eV 
[3T][52]; the resulting relative positions of the first mo- 
ments and the large conduction band overlap indicate 
again a good applicability of the CPA. In contrast to the 
previous material systems, the spin-orbit coupling is one 
order of magnitude smaller in the nitride compounds and 
does not influence the results for the bowing. 

The E g (x) results for Ca^Ali^N are depicted in 
Fig.[l0| The numerical parameters for the CPA and the 
supercell calculations were chosen identical to those of 
the previous sections. In all three models, we can clearly 
identify a crossover in the bowing between x — 0.3 and 
x = 0.2. In this region, the character of the band gap 
changes from a behaviour which is strongly influenced by 
the indirect A1N material to a direct band gap behaviour 
dominated by GaN. 

While the CPA results coincide very well with the su- 
percell results for the GaN-dominated side at large x, 
they overestimate the bowing on the Al-rich side. In con- 
trast, the VCA, which shows a piecewise linear behaviour 
with two different slopes, can reproduce the supercell re- 
sults for high Al contents quite well, but deviates for 
0.7 > x > 0.3. It additionally should be noted that the 
CPA results in the Al-rich part have larger error bars, 
as the determination of the band gap is afflicted with a 
larger uncertainty in this range due to a larger broaden- 
ing of the corresponding spectral functions (not shown). 

In [22] , Vurgaftman et al. report bowing parameters of 
b r = 0.05-0.53 eV for the T-valley of cubic Ga^Ali^N 
from theory. By additionally taking several experimental 
results into account (which obviously render larger bow- 
ing parameters), they recommend an approximate value 
of br ~ 0.7 eV. By only fitting the T-valley bowing, i. e. 
only taking the values for x > 0.3 into account, aug- 
mented by the A1N energy difference at T as boundary 
value, we obtain bowing parameters of br = (0.37 ±0.02) 
eV from the supercell calculations and br = (0.4±0.1) eV 
from the CPA. Consequently, the CPA and supercell val- 
ues agree within the error range and our results are in rea- 
sonably good agreement with the literature values. More 
recent results from DFT+LDA calculations [35] yield a 
value of br ~ 0.5 eV and predict the crossover at about 
x = 0.4. However, it should be noted that these results 
suffer from the usual underestimation of band gaps in 
DFT+LDA. The pure GaN and A1N band gaps in their 
calculations are obtained as 1.93 eV and 3.23 eV respec- 
tively. Consequently, they strongly deviate from our in- 
put values of 3.26 eV for GaN and 5.346 eV for A1N 
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FIG. 7: (Color online) Band structure (left) and DOS and first moments (right) of InAs and GaAs, calculated in the sp 3 
EBOM. Again, the energies of InAs are shifted by a valence band offset A_E v b = 0.85 eV with respect to GaAs. See appendix 
for details. 
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FIG. 8: (Color online) Concentration-dependent band gap 
E s (x) of InzGai-zAs including spin-orbit coupling, calcu- 
lated in the VCA (blue), Wannier CPA (red) and supercell 
method (black). 

(see appendix), as our CPA/supercell + EBOM model 
allows for the usage of arbitrarily exact boundary values 
at x — and x = 1. 

For the sake of comparison, we summed up the results 
for the bowing parameters in Tab.||J As already stated, 
the given error ranges only account for the reading accu- 
racy, finite size effects and non-parabolicity of E g (x) but 
cannot reflect the reliability of the input band structure 
parameters for the pure materials. 



IV. CONCLUSION AND OUTLOOK 

In this paper, we showed that the electronic properties 
of substitutional semiconductor alloys with an underly- 
ing zincblende structure of the type A^Bi-^C can well 



TABLE I: Summary of the results for the bowing parameter 
b in eV. In case of Ga^Ah-^N, the value solely refers to the 
bowing of the F-valley. 



Material 


Supercells 


Wannier CPA Literature 


Cd^Zm-^Se 
In^Gai-^As 
Ga^Ali-^N 


0.71 ±0.09 
0.39 ±0.04 
0.37 ±0.02 


1.0 ±0.1 0-1.26 
0.49 ± 0.05 0.32-0.46 
0.4 ±0.1 0.05-0.7 



be described with empirical TB models. We presented 
a combined theoretical approach, starting from a multi- 
band TB model with a realistic dispersion and bandwidth 
throughout the whole Brillouin zone. The density of 
states and the band gap of the disordered system can 
then either be determined by exact diagonalization of 
a large supercell with a large number N of microscopi- 
cally distinct configurations, or by combination of the TB 
model with the coherent potential approximation (CPA) 
and/or the virtual crystal approximation (VCA). 

Using the supercell results for Cd^Zni-^Se as refer- 
ence, we gave a careful quantitative and qualitative anal- 
ysis of the scope of validity of the CPA and VCA, espe- 
cially with regards to the calculation of the concentration 
dependent band gap E e (x) of the alloy. While the VCA 
failed over the whole concentration range, the CPA also 
turned out to be not accurate enough for the Cd^Zni-^Se 
system under consideration, although the proper choice 
of the basis set could significantly reduce the discrepancy. 

We then applied our TB model to two further different 
alloy systems, namely the III-V alloy In^Gai-^As and 
the Ill-nitride system Ga^Ali^^N. For both systems, the 
CPA gave good results. In case of In^Gai-^As, the CPA 
and the supercell calculations yielded bowing parame- 
ters in good agreement with literature values from ex- 
periments. For Ga^Ali-^N the band gap bowing showed 
a crossover behaviour between x = 0.3 and x = 0.2, due 
to the fact that cubic GaN has a direct energy gap at the 
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FIG. 10: (Color online) Concentration-dependent band gap 
E g (x) of Ga^Ali-^N, calculated in the VCA (blue), Wannier 
CPA (red) and supercell method (black). 



combined with the ETBM and especially the EBOM, the 
application to further material systems (with disorder ei- 
ther in the cations or in the anions) will be an interesting 
task for the future, ideally alongside actual experimen- 
tal data. Furthermore, the same calculation scheme can 
be transferred to alloy systems with underlying wurtzite 
structure by using a suitable TB Hamiltonian [15j for 
direct as well as indirect band gap materials. 

Our supercell approach is also applicable to disordered 
low-dimensional structures, as shown in [I], [IB]. There- 
fore, in principle also disordered nanowires or superlat- 
tices can be investigated. As the CPA is exact only in the 
limit of infinite dimensions, it may turn out that the CPA 
+ EBOM is less reliable when applied to low-dimensional 
systems, so that suitable extensions of the CPA must be 
used. 



BZ center, while the conduction band minimum of cubic 
A1N is located at X. This crossover was reproduced in 
all three models. The CPA and the supercell calculations 
were able to reproduce the T-valley bowing in satisfac- 
tory agreement with the literature. On the Al-rich side 
(x < 0.3), the CPA understimated the band gap when 
compared to the supercell approach, while the VCA was 
in surprisingly good agreement. 

For the sake of completeness, it should be emphasized 
that the computational costs of the CPA are far smaller 
than in the supercell case. If the number of bands has 
to be augmented or properties far from the band edges 
become relevant, the supercell approach can quickly be- 
come infeasible, as the calculation time scales with the 
cube of the dimension of the Hamiltonian matrix. 

As the supercell calculations and, under certain condi- 
tions outlined in this paper, also the CPA can give good 
results for the concentration-dependent band gap when 



Appendix 

The following tables give the material parameters used 
throughout this paper, including the sources from the 
literature. 

The band structures in the sp 3 EBOM can be fitted to 
the conventional lattice constant a, the spin-orbit split- 
ting A so , the effective conduction band mass m c , the 
Luttinger parameters 71,72,73, and to a set of energies 
at the T-point and the A-point (denoted by the usual 
single group notation), where the energy gap E g is given 
as — in case of a direct band gap. Additionally, 
a valence band offset AE v t, is incorporated in order to 
account for the relative energetic position of the two con- 
stituents. 

These parameters are uniquely connected to the non- 
vanishing TB matrix elements from Eq. ([I]) of the present 
paper by the equations (7)-(17), (33) and (45)-(53) of 
Ref.H 



Appendix A: Material parameters CdSe and ZnSe 



Parameter 


CdSe 


ZnSe 


a 


(A) 


6.078 


33 


5.668 




A so 


(eV) 


0.41 


03] 


0.43 


33] 


m c 


(mo) 


0.12 


03] 


0.147 


34j 


T^C T^V 

1 1 1 15 


(eV) 


1.76 


03 


2.82 


33J 


xi 


(eV) 


2.94 


[36] 


4.41 


36J 


XI 


(eV) 


-1.98 


03 


-2.08 


03 


xi 


(eV) 


-4.28 


03 


-5.03 


03 


7i 




3.33 


02 


2.45 


02 


72 




1.11 


02] 


0.61 


34] 


73 




1.45 


03] 


1.11 


34j 


A£ vb 


(eV) 


0.22 


03 





03] 
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Appendix B: Material parameters In As and GaAs 



Parameter 


InAs 


GaAs 




(A) 


6.058 [37] 


5 653 




A 

^so 


\ e v ) 


0.39 




U.04: 




m c 


(Tfln) 


0.022 




0.067 


I2TI 


1 1 1 15 


(eV) 


0.417 


m 


1.519 


02 


X{ 


(eV) 


2.28 


03 


2.18 


38. 


xi 


(eV) 


-2.42 


03 


-2.8 


38J 


xi 


(eV) 


-6.64 


m\ 


-6.7 


38J 


7i 




20.0 


m 


6.98 


21J 


72 




8.5 


m\ 


2.06 


02 


73 




9.2 


m\ 


2.93 


ED 


AE vh 


(eV) 


0.85 


[27] 








Appendix C: Material parameters GaN and A1N 



Parameter 


GaN 


A1N 


a 


(A) 


4.50 


m 


4.38 


31 


A so 


(eV) 


0.017 


m\ 


0.019 


0H 


m c 


{mo) 


0.15 


0B 


0.25 


02 


1 1 1 15 


(eV) 


3.26 


0B 


5.84 


30 


xi 


(eV) 


4.43 


00] 


5.346 


30 


xi 


(eV) 


-2.46 


00] 


-2.315 


30 


xi 


(eV) 


-6.30 


m 


-5.388 


30 


71 




2.67 


0B 


1.92 


02 


72 




0.75 


m\ 


0.47 


02 


73 




1.10 


0B 


0.8* 




A£ vb 


(eV) 


0.8 


M 





02 



*This parameter has been adjusted by hand, as the orig- 
inal value of 73 =0.85 used in |31] leads to an erroneous 
curvature in the Y—K direction. 



* Electronic address: dmourad@itp.uni-bremen.de 
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